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Abstract 

P-Stars are a new class of compact stars made of up and down quarks in /^-equilibrium with 
electrons in an Abelian chromomagnetic condensate. We show that P-Stars are able to account 
for compact stars with R < 6 Km, as well as stars with radius comparable with canonical Neutron 
Stars. We find that cooling curves of P-Stars compare rather well with observational data. We 
suggest that P-Matter produced at the primordial deconfinement transition is a viable candidate for 
baryonic Cold Dark Matter. Finally, we show that P-Stars are able to overcome the gravitational 
collapse even for mass much greater than 10 6 Mq. 
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One of the most intriguing aspect of modern astrophysics is the possible existence of 
Strange Stars which are made entirely of deconfmed u, d, s quarks. The possible 

existence of Strange Stars is a direct consequence of the so called Strange Matter hypothesis, 
namely that u, d, s quarks in equilibrium with respect to weak interactions could be the 
true ground state of hadronic matter 

Recently, numerical lattice results in QCD 0, |(| suggested that the gauge system gets 
deconfmed in strong enough Abelian chromomagnetic field Af.(x) = 5 a ' 3 8k,2 %i H. This leads 
us to consider a new class of compact quark stars made of almost massless deconfmed up 
and down quarks immersed in a chromomagnetic field in /3-equilibrium. As we will discuss 
below, it turns out that these compact stars, which will be referred to as P-Stars, are more 
stable than Neutron Stars whatever the value of y/gH, g being the color gauge coupling. 
To investigate the structure of P-Stars we needs the equation of state appropriate for the 
description of deconfmed quarks and gluons in an Abelian chromomagnetic field. The quark 
chemical potentials turn out to be smaller that the strength of the chromomagnetic field 
yfgH. Thus, we see that up and down quarks are in the lowest Landau level with energy 
£o,p3 = \ps\, where p is the quark momentum. It is now easy to find the thermodynamical 
potential (at T = 0). We find for the energy density: 



4tt 2 ^ ^ u ^ a> An 2 32tt 2 
while the pressure P is given by: 



9 H (fl+fd) + TZi + . (1) 



where /// ( / = d,u,e) denote the chemical potential. Moreover, we have the following 
constraints: 

A* e + /J-u = P — equilibrium (3) 



-^ n u — -^ n d = n e charge neutrality (4) 



where: 

1 „ _ 1 „ _ /C 



n„ 



gHu u , n d = —gHu d , n e = — e - . (5) 



We obtain the following Equation of State: 
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It is interesting to observe that the sound speed: 

2 _ dP 1 m 

is quite close to the causal limit vs = 1, for /Z e = ^ e /^gH < 1, leading to a very stiff 
equation of state. 

To study some properties of a star we consider it to be a spherical symmetric object, cor- 
responding to a non rotating star. The stability of the star is governed by the general- 
relativistic equation of hydrostatic equilibrium for a spherical configuration of quark matter, 
which is the Tolman-Oppenheimer-Volkov equation p|: 

dP GM(r)e(r) u , 4 7 rr 3 P(r) 1 1 + 



1 + iift„\ ] : 2GMM ' ( 8 ) 



dr " r 2 1 M(r) J i _ 2GM ^ ' 

dM 



Anr 2 e(r). (9) 
dr 

These equations can be solved numerically for a given central density e c and obtain e(r) and 
P(r). The radius of the star is determined by : 

P(r = R) = , (10) 

and the total mass by: 

M = M(r = R) . (11) 

In Figure 1 we display the gravitational mass M versus the star radius R for different values 
of y/gH. At fixed value of the chromomagnetic condensate the mass and radius of the star 
can be thought of as function of the central energy density e c . We see that, as in Strange 
Stars, the mass first increases with e c until reaches a maximum value where dM }- £c " > = 0. 
Increasing further e c leads in a region where dj ^ £c - > < and the system becomes unstable. 
Moreover, we see that there is no lower limit for the radius R. Indeed, for small mass we 
find: 

R ~ M5 . (12) 

Interestingly enough, Fig. 1 shows that there are stable P-Stars with M ^ Mq and R ^ 
6 Km. This region of the M — R plane could be relevant for the recently observed compact 
star RXJ1856. 5-3754 3,3- In fact > RXJ1856. 5-3754 has been observed with Chandra and 
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FIG. 1: Gravitational mass M plotted versus stellar radius R for P-stars for different values of 
\JgH. Thick solid line corresponds to M-R curve for RXJ1856. 5-3754 obtained solving Eq. (|14jl 
with R°° = AAKm and R>R 1 = 3GM . 

XMM-Newton [10], showing that the X-ray spectrum is accurately fitted by a blackbody 
law. Assuming that the X-ray thermal emission is due to the surface of the star, the authors 
of Ref. [nj] found for the effective radius and surface temperature: 

R°° ~ 4.4 -^—Km , T°° ~ 63 el/, (13) 
120pc ' ' v ; 

where: 



R-= , R ,i~ = rJi-ggg. (14) 

l-y 2GM V R 



R 

Assuming R°° = 4 .4 Km (i.e. d = 120 pc) we can solve Eq. ()14j) for the true radius R. We 
also impose that R > R 7 = 3GM, for it is well known that the last circular photon orbit 
for the Schwarzschild geometry is at Rj. Remarkably, Fig. 1 shows that there are stable 
P-Star configurations which agree with observational data. However, it should be stressed 
that in the observed spectrum there is an optical emission in excess over the extrapolated 
X-ray blackbody. Indeed, optical data are also fitted with a blackbody model yielding an 
effective radius > 17 Km 12 p |lo||- We feel, however, that the seven data points in 
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the optical range could be interpreted differently. For instance, the optical data could be 
interpreted as syncroton radiation emitted by electrons with spectrum: 

N(E) = kE~ v , 2. 10~ 3 KeV < E < 7. KT 3 KeV. (15) 

We find a quite good fit with r] ~ 0.71. If this is the case the radiation in the optical range 
should display a polarization II = — 0.56. 

Beside, there are stable P-Stars with mass and radius comparable to those of Neutron Stars. 
In this case, we show that P-Stars are more stable than Neutron Stars. To see this, we need 
the baryon number: 

A = An f R ^ = r * dr ( 16 ) 

' j _ 2GM(r) 



where n b = \{n u + n^). Thus, the binding energy per nucleon is : 



B = m n -^, (17) 



m n being the nucleon mass. We find for stars with limiting mass: 



B ~ 0.125 GeV , A ~ 1.5 10 57 , — — ~ 1.1 

Mq 

B ~ 0.386 GeV , A ~ 2.9 1 57 , — ~ 1.4 

Mq 



(18) 



These results should be compared with the binding energy per nucleon in the case of limiting 
rn 

mass Neutron Star [7], B^s — 100 MeV, showing that P-Stars are gravitationally more 
bounded than Neutron Stars. 

Moreover, it turns out that P-stars are also more stable than strange stars . Indeed, 
it is known that there is extra energy of about 20 MeV per baryon available from the 
conversion of matter to strange matter So that, even the binding energy per nucleon 
for strange stars, B$s — 120 MeV, is smaller than the value in Eq. (fTHj) . This shows 
that gravitational effects make two flavour quark matter in the chromomagnetic condensate 
globally more stable than nuclear matter or three flavour quark matter. However, one could 
wonder if locally, let's say in a small volume of linear size a ~ ^-/fn n ~ 1 fermi, two flavor 
quark matter immersed in the chromomagnetic condensate could convert into neutrons or 
deconfined quarks of relevance for strange stars. In order to see how far this conversion 
will proceed, we need to estimate the appropriate transition rates. Let us first consider the 
decay into nuclear matter. In this case we must rearrange three quarks, which are in the 



lowest Landau levels, into configurations of quarks in nucleons. In doing this, we get a first 
suppression factor a due to color mismatch. In fact, the quarks into the lowest Landau 
levels are colorless due to presence of the chromomagnetic condensate. This suppression 
factor cannot be easily evaluated. However, in general we have that a < 1- Another 
suppression factor arises from the mismatch of the quark wavefunctions in the directions 
transverse to the chromomagnetic field. This results in a factor (^=2=^ for each quarks. 
Finally, to arrange three quarks into a nucleon we need to flip one quark spin, which increases 
the energy by a factor y/gH at least. So that we have for the transition rate: 

p ~ a \vm) e • (19) 

where T is the star core temperature. For a newly born compact star the core temperature 
is of order of several MeV and it is rapidly decreasing. On the other hand, the typical values 
of y/gH is of order of several hundred of MeV. As a consequence we have: 

P < 1.5 10" 46 , (20) 

which implies that the decays of two flavor quark matter immersed in the chromomagnetic 
condensate into nuclear matter is practically never realized. In addition, the conversion into 
three-flavor quarks is even more suppressed. Indeed, first we need to flip a large number 
of quark spin to fill the Fermi sphere, then the quarks at the top of the Fermi sphere are 
allowed to decay into strange quarks. We conclude, thus, that P-stars, once formed, are 
absolutely stable. 

The binding energy is the energy released when the core of an evolved massive star 
collapses. Actually, only about one percent of the energy appears as kinetic energy in the 
supernova explosion jll|. Now, if we admit that supernova explosions give rise to P-Stars 
instead of Neutron Stars, then we may completely solve the supernova explosion problems, 
for there is an extra gain in kinetic energy of about 1 — 10 foe (1 foe = 10 51 erg). Then we 
could eventually suppose that Pulsars are P-Stars and not Neutron Stars. Further support 
to this point of view comes from cooling properties of P-Stars. 

Let us assume, for simplicity, stars of uniform density and isothermal. The cooling equation 
is: 

C v ^ = -(L u + L 7 ) , (21) 



6 



where L v is the neutrino luminosity, L 1 is the photon luminosity and Cy is the specific heat. 
Assuming blackbody photon emission from the surface at an effective surface temperature 
T s we get: 

L 7 = AnR 2 a SB Tl, (22) 

where ctsb is the Stefan — Boltzmann constant. Following Ref. we assume that the 
surface and interior temperature are related by: 

T 

Y = 10" 2 a , 0.1 < a < 1.0 . (23) 

The dominant cooling processes by neutrino emission are the direct /3-decay quark reac- 
tions 3,Q : 

d — ► u + e + V e , , u + e — > d + u e , . (24) 

Let \M\ 2 be the squared invariant amplitude averaged over initial spin and summed over 
final spins. A standard calculation gives: 

|M| 2 = Q4G 2 F cos 2 6 c (p-q')(p'-q) , (25) 

where the quark momenta p and p' are directed along the chromomagnetic field direction, q 
and q' are the electron and neutrino momentum respectively . Now energy and momentum 
conservation give: 

p = (E p ,p 3 ) , \p 3 \ ~ fi d ; p' = (Ep^pz) , \p' 3 \ ~ /i u 
q = (E g ,q 3 ,q ± ) , \q 3 \ ~ /i e , |g ± | ~ T (26) 
q' = (E g ,,c7), \q\~T. 

So that we find: 

\M\ 2 ~ 32G 2 F cos 2 9 c E p E p , E q E g , (— ) 2 . (27) 

Using standard methods for calculating phase space integrals in degenerate Fermi sys- 
tems [15], we find the following neutrino emissivity: 

1371 Fir 



30240 V 2 
Thus, the neutrino luminosity is: 



G F cos 2 9 C (7i e ) 2 ^gH T 8 . (28) 



3 ' 18 10 V Tg Mo 7 TGeV ' (29) 
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FIG. 2: P-star cooling curves for M = 0.6 Mq and M = 2.7 Mq and three different values of the 
parameter a in Eq (|2«-i|) . Experimental data for the effective surface temperature have been taken 



from Refs. 
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16, 17]. 



where T 9 is the temperature in units of 10 °K, and Eq = 2.51 10 \ gr/ cm is the nuclear 



density. Finally, the specific heat is given by: 



0.92 10 55 T ( M £ ° 





Mq e K lGeV 



(30) 



We obtain the thermal history of a P-Star by integrating Eq. (}2~T]) with T 9 



1.4. 



In Figure 2 we report our cooling curves for three different values of the parameter a in 
Eq. ()23)1 . It is worthwhile to stress that the cooling curves are almost independent on the 
P-Star mass. Moreover, there is a weak dependence on a up to age r ~ 10 3 years. We 
compare our theoretical cooling curves with Pulsar data taken from Ref . Ilfil . We also report 
RXJ 1856. 5-3754 effective surface temperature assuming r ~ 10 5 years [lOj, and the surface 



effective temperature upper limit for supernova remnant 3C58 |17j. We stress that the upper 
limit has been obtained with R°° = 12 Km. Assuming an effective radius smaller than the 
Neutron Star canonical value should increase the above upper limit. In any case we see that 
the agreement between theoretical cooling curves and observational data is quite satisfying. 
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Long time ago Witten [4] conjectured that quark nuggets could be copiously produced 
during the cosmological deconfmement phase transition, providing an explanation for the 
Dark Matter in the Universe. Whether or not Strange Matter j3| is absolutely stable, the 
authors of Ref. Q argued that a .urge lump of Strange Matter evaporates on.y its outer 
layer, so that lumps large enough could survive until today. In fact, a detailed analysis [19J 
showed that minimum baryon number of lump that survives is roughly 10 52 , restricting 
severely the viability of Strange Matter as a candidate of Dark Matter. On the other hand, 
we now show that this limitation is not shared by P-Matter, i.e. lumps of up and down 
quarks in equilibrium with respect to the weak interactions immersed in a chromomagnetic 
condensate. 

Recently, it has been suggested that strong magnetic fields of order 10 19 G, corresponding to 
y/gH ~ 1 GeV, are naturally associated with the QCD scale j^j]. Moreover, it is believed 
that large magnetic fields might be generated during cosmological phase transitions [21]. So 
that it is conceivable that P-Matter with chromomagnetic condensate y/gH ~ 1 GeV could 
have been copiously produced at the cosmological deconfmement transition. The physical 
properties of P-Matter are very similar to those of Strange Matter. Let A indicate the 
baryon number of a lump of P-Matter, then we find: 

M ~ (2.5GeV) A,R^ , (31) 

where ~ (248 MeV) 3 . Now, the baryon number of P-Matter can change only through the 
emission or absorption of neutrons. Moreover, we estimate the binding energy / > \fg~H- 
As for Strange Matter, the neutron cross section is almost geometric: 

a ~ 4nR 2 f n ~ f n avA 2/ \ (32) 

where cr ~ 0.79 lO^MeU" 2 and /„ < ^ 0.04. Following Ref. Q we find: 

^ = - r(A) , r(A) * ^ 2 T 2 e-i f n a A 2 / 3 . (33) 

The evaporation rate of a lump of P-Matter when the Universe had a temperature Tu is 
given by: 

dA dA dt dA I 45" 

(34) 



dT v dt dTu dt V 172tt 3 G T 3 ' 
Integrating Eq. (JB1J), assuming the conservative value /„ ~ 0.1, with typical y/gH ~ 1.0 GeV 

and initial temperature Tu ~ T c ~ 150 MeV, we find that the minimum size lump of P- 

Matter that survives until today has a baryon number A min ~ 10 42 . This number should be 
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FIG. 3: P-star gravitational mass M versus stellar radius R for \/gH = 10~ 4 GeV. Dash-dotted 
line corresponds to the Schwarzschild radius i?g = 2GM . 

compared with A c , the mean baryon number per horizon at T c assuming that the Universe is 
closed by baryons: A c ~ 10 48 . We see that A min is much smaller than A c , making P-Matter 
a viable candidate for Cold Dark Matter. 

The most interesting and intriguing aspect of our theory is that P-Stars are able to overcome 
the gravitational collapse even for mass much greater 10 6 Mq. To see this, we note that on 
dimensional ground we have: 

1 » 1 



M 



f(Sc 



R 



G^gH Jy ^> " G 1 / 2 gH 
where e c = e c j (gH) 2 . As a consequence we get: 



(35) 



2GM 



f(e c 



h(e c 



(36) 



From Equation (|35jl we see that by decreasing the strength of the chromomagnetic conden- 
sate we increase the mass and radius of the star. However, the ratio 2 R M depends only on e c . 
It turns out that the function h(x) defined in Eq. ()36p is less than 1 for any allowed values of 
e c . Indeed, in Fig. 3 we shows stable configurations with M ~ 10 7 Mq and R > 1.3 Rg, Rs 
being the Schwarzschild radius. Thus, we infer that our P-Stars do not admit the existence 



10 



of an upper limit to the mass of a completely degenerate configuration. In other words, our 
peculiar equation of state of degenerate up and down quarks in a chromomagnetic conden- 
sate allows the existence of finite equilibrium states for stars of arbitrary mass. 
Let us conclude by stressing the main results of this paper. We found that P-Stars are able 
to describe very compact stars as well as stars in the Neutron Star region. We showed that 
the cooling curves compare reasonable well with observational data. Moreover, we found 
that P-Matter could be a candidate for baryonic Cold Dark Matter. However, the most 
striking result is that P-Stars are able to halt the gravitational collapse. So that we may 
conclude that P-Stars are challenging Neutron Stars and Black Holes. 
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